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Abstract 

Stimulated by a recent development of the universal seesaw mass matrix 
model, the evolutions of the gauge and Yukawa coupling constants are inves- 
tigated under the gauge symmetries SU(3) c xSU(2)£xSU(2)#xU(l)y. Espe- 
cially, an investigation is made as to whether this evolution can constrain the 

necessary intermediate scales in these types of models and its viability. 
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I. INTRODUCTION 



Recently, the so-called "universal seesaw mass matrix model" [ffj has been revived as 
a model which gives a unified description of masses and mixings of the quarks and leptons. 
The "seesaw mechanism" was first proposed j|] in order to answer the question of why 
neutrino masses are so invisibly small. Then, in order to understand that the observed quark 
and lepton masses are considerably smaller than the electroweak scale = (<f)° L ) = 174 GeV, 
the mechanism was applied to the quarks ffl. However, the observation of the top quark of 
1994 P| aroused a doubt on the validity of the seesaw mechanism for the quarks because the 
observed fact m t ~ means that M^ttlr is of the order of one in the seesaw expression 
Mf ~ m^Mp^mR. On the contrary, it has recently been found |2||| that the model can give 
an interpretation for the question of why only top quark acquires a mass of the order of A^ 
if we take an additional condition detM^ = for up-quark sector. 

In the universal seesaw mass matrix model, the mass matrix for fermions (/, F) is given 

by 



M 



m 



L 



I m L \ ( Z- » 



i kZr \Yp j 



where /j (fermion sector names / = u,d,i>,e; family numbers % = 1,2,3) denote quarks 
and leptons, F% denote hypothetical heavy fermions F = U, D, N and E correspondingly to 
/ = u,d,v and e, and they belong to f L = (2, 1), f R = (1, 2), F L — (1, 1) and F R = (1, 1) 
of SU(2)lxSU(2)^. The matrices Zl, Zr and Yp are those of the order of one. The 3x3 
matrices (~ itlq = Al) and tur (~ Km® = Ar) are symmetry breaking mass terms of 
SU(2)i and SU(2)#, respectively, and those have common structures independently of the 
fermion sector names /. Only Mp (~ Xtuq = As) has a structure dependent on the sector 
name /. For the case A ^> k ^> 1, the mass matrix (1.1) leads to the well-known seesaw 
expression 

M f ~ m L Mp l m R . (1.2) 
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In contrast to the case (1.2), for the case with the additional condition 



detM F = , (1.3) 

on the up-quark sector (F = U), one of the heavy fermions F, (say, F 3 ) cannot acquire a 
mass of the order of A5 = Xmo, so that the seesaw mechanism does not work for the third 
fermion. Therefore, the mass generation at each energy scale is as follows: First, at the 
energy scale /i = As, the heavy fermions F, except for U3, acquire the masses of the order of 
As- Second, at the energy scale \x = An, the SU(2)r symmetry is broken, and the fermion 
ur3 generates a mass term of the order of Ar by pairing with Ul3- Finally, at \i = Al, the 
SU(2)i symmetry is broken, and the fermion Ui3 generates a mass term of the order by 
pairing with Urz- The other fermions / acquire the well-known seesaw masses (1.2). The 
scenario is summarized in Table |. We regard the fermion pair (ul 3 , Ur 3 ) as the top-quark 
state. Thus, we can understand why only top quark t acquires the mass m t ~ O(m^) @,[|. 

On the other hand, for the neutrino mass generation, at present, we have the following 
two scenarios as summarized in Table |TJ. One (Scenario A) is a trivial extension of the 
present model: we introduce a further large energy scale A uS in addition to As, and we 
assume that M F ~ A s (F = U,D,E), while M N ~ A uS (A uS A s ). Another scenario 
(Scenario B) || is one without introducing such an additional energy scale. The neutral 
heavy leptons are singlets of SU(2)ixSU(2)^ and they do not have U(l)-charge. Therefore, 
it is likely that they acquire Majorana masses Mm together with the Dirac masses Mr, = 
at fi = As- For example, we assume Mm = Md [0. Then, the neutrino mass matrix for the 
conventional light neutrinos is given by M v = m^Mjfm^, so that the masses m v are given 
with the order of 

A| 1 A L A R 

~A S ^ ~k~As~ ' (L4) 

In order to explain the smallness of m v , the model requires that the scale Ar must be 
extremely larger than A^ (for example, k = A^/A^ ~ 10 9 0). This scenario seems to be 
very attractive from the theoretical point of view, because we can explain the mass hierarchy 



of the quarks and leptons by the three energy scales A^, Ar and As only. On the other 
hand, in the scenario A, there is no constraint on the value of n (however, the value must 
be larger than ~ 10 because of no observation of the right-handed weak bosons Wr at 
present), so that the model allows a case with a lower value of Ar. Since we can expect 
abundant new physics effects for the case of k ~ 10 ||, the case is also attractive from the 
phenomenological point of view. 

One of the purposes of the present paper is to see whether a study of the evolutions of the 
gauge coupling constants of SU(3) c xSU(2)LxSU(2) j RxU(l)y and of the Yukawa coupling 
constants in the universal seesaw mass matrix model can give any hint on the value of the 
intermediate energy scale A R or not. For example, Shafi and Wetterich || and Rajpoot 
T0[1 have considered an O(10) model and an SO(10) model, respectively, with the symmetry 



breakings SO(10)-> SU(3) c xSU(2) L xSU(2) i? xU(l) y at /i = A GUT and SU(2) fi ->U(l) fl at 
/i = Ar, and they have demonstrated that the model with Aqut ~ 10 19 GeV and Ar ~ 10 9 
GeV is consistent with the low energy phenomenology. The value A^ ~ 10 9 GeV is favorable 
to the scenario B for neutrino masses. However, in the present model, since there are many 
new fermions F above the intermediate energy scale As, their conclusion cannot be applied 
to the present seesaw mass matrix model straightforwardly. 

On the other hand, a phenomenological study of the universal seesaw mass matrix model 
for the quark masses and the Cabibbo-Kobayashi-Maskawa (CKM) Jll| matrix parameters 
has successfully been given by Fusaoka and the author ||. In order to give explicit numerical 
predictions, they have used some working hypotheses that I will use here as well. 

(i) The matrices Zl and Zr, which are universal for quarks and leptons, have the same 
structure: 

Z L = Z R = Z = diag(zi, z 2 , z 3 ) , (1.5) 

with z\ + z\ + z\ = 1 , where, for convenience, we have taken a basis on which the matrix Z 
is diagonal. 

(ii) The matrices Yp, which have structures dependent on the fermion sector f = u,d, v, e, 



take a simple form [(unit matrix) + (a rank one matrix)]: 



Y f = 1 + 3b f X 



(iii) The rank one matrix X is given by a democratic form 



(1.6) 



'lll^ 



1 1 1 
1 1 1 



(1.7) 



on the family-basis where the matrix Z is diagonal. 

(iv) In order to fix the parameters z iy we tentatively take b e = for the charged lepton 
sector, so that the parameters Z{ are given by 

1 



Zl 



Z2 



Z 3 



By taking b u = —1/3 (then detMjj = 0), they have obtained the following top-quark mass 
enhancement without the suppression factor k/X 

1 



(1.9) 



together with the successful relation m u /m c ~ 3m e /4m^. Furthermore, by taking bd = 
—e' l/3d (f3 d = 18°), they have succeeded in giving the reasonable values of the CKM matrix 
parameters together with the reasonable values of the quark mass ratios (not only mf/rrij, 
mf/rrij, but also mf/rrij) with keeping the value of the parameter (m /t/A)/ in (m /t/A) M = 
(mon/X)d- However, in order to fit the quark mass values (not the ratios) to the observed 
quark mass values at /i — m z , they have taken the parameter (m /t/A)/ as 

(m K/X) u \ 



R{m z ) = 



(m K/X) ( 



~ 3 



;i.io) 



H=mz 



It seems to be natural to consider that all Yukawa coupling constants become equal between 
quarks and leptons at a large energy scale A Y u- Therefore, another one of the purposes of 
the present paper is to see whether such a factor 3 can be understood by the difference of the 
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evolutions of the Yukawa coupling constants between quarks and leptons from the energy 
scale /i = A Y u to p = rriz- 

In Sec. II and Sec. Ill, we investigate evolution of the gauge and Yukawa coupling con- 
stants, respectively, under the gauge symmetries SU(3) c xSU(2)LxSU(2)pxU(l)y at one- 
loop. We will conclude that it is possible to find the energy scale Ayu at which R{n) takes 
R = 1 only for a model with a value k < 10 2 . Although in Sec. II and Sec. Ill we consider 
the case that the symmetries SU(3) c xSU(2)£,xSU(2)pxU(l)y are unbroken for the region 
> As, in Sec. IV, we investigate a case that the symmetries SU(3) c xU(l)y are embed- 
ded into the Pati-Salam symmetry [12| SU(4)ps at fi > As, so that we consider the case 
of SU(4) P5 xSU(2) L xSU(2) fi in the region A s < pi < A GUT . We will find that the model 
predicts A R ~ 5 x 10 12 GeV, A 5 ~ 3 x 10 14 GeV, and A GUT ~ 6 x 10 17 GeV. Finally, Sec. V is 
devoted to the conclusions and remarks. We will find that there is no model which satisfies 
Ayu = Aqut- 

II. EVOLUTION OF THE GAUGE COUPLING CONSTANTS 

The gauge symmetries SU(3) c xSU(2) i xSU(2)pxU(l)y are broken into the gauge sym- 
metries SU(3) c xSU(2) i xU(l)y/ at ji = A R . The electric charge operator Q 

Q = lf + I* + ±Y, (2.1) 

at fi > An is changed into 

q = ^ + \y , (2.2) 

in the region Al < jj, < Ar. Hereafter, we call the regions Al < ft < Ap, Ar < fi < As and 
A5 < p < A x (A x = Aqut or A x = A Y u) Regions I, II and III, respectively. 

The evolutions of the gauge coupling constants gi at one-loop are given by the equations 

4«i(^) = -^-6<o<(^) > (2-3) 
at 2n 
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where 014 = gf/4n, t = ln/i and the coefficients bi are given in Table [TV]. (Note that the heavy 
fermions F L and F R except for Ul3 and U R3 are decoupled for /1 < A s and the fermions u R3 
and ^3 are decoupled for /i < A R .) The boundary conditions at /i = Al and /i = A R are 
as follows: 

^(Al) = « l 1 (Al) + ^T 1 (Al), (2.4) 

and 

l a 'r\A R ) = a R \A R ) + ^\A R ) , (2.5) 

respectively, correspondingly to (2.2) and (2.1), where the normalizations of the U(l)y 
and U(l)y gauge coupling constants have been taken as they satisfy a[ = = a 3 in the 
SU(5) grand-unification limit and a.\ = = a R = a 3 in the SO(10) grand-unification 
limit, respectively. For convenience, we use the initial values at /x = rriz instead of those 
at /J, — A L in the region I: a[(mz) = 0.01683, a^mz) = 0.03349 and a 3 (m z ) = 0.118. 
The values of a[ and ai have been derived from ]13[ a em {mz) = (128.89 ± 0.09) -1 and 
sin 2 9w = 0.23165 ± 0.00024. The value of a 3 has been quoted from Ref. |b||. We illustrate 
a typical case with A R = 10 5 GeV and a R (A R ) = l/4ir in Fig. 1. 

In the numerical study, we have taken the value of the parameter A R /As = k/X as 

k/X = A R /A S = 0.02 , (2.6) 

which has been obtained form the observed value of the ratio m c /mt in Ref. Q. Although 
the value (2.6) has been obtained on the model with the specific matrix forms (1.5) - (1.7), 
the order of the value (2.6) will be valid for any other seesaw model with detMu = because 
in such a model the value of k/X is given by the order of mb/m t . 

As seen in Fig. 1, the U(l) coupling constant ai(/x) becomes rapidly strong in the region 
III (/i > As) because the heavy fermions F become massless in the region III. We consider 
that the unification energy scale Ayu of the Yukawa coupling constants must be lower than 
an energy scale Af 3 at which Q!i(/x) becomes infinity. This condition will impose a strong 



restriction on the possible Ay;y-search as we discuss in the next section. Of course, in the 
grand unification scenario, the U(l) symmetry will be embedded into a grand unification 
symmetry G before the U(l) coupling constant bursts. Such a case will be discussed in 
Sec. IV. 



III. EVOLUTION OF y L y R /y s 

The 3x3 matrices tjil, rn R and Mp are given in terms of the vacuums expectation values 
vl = V2((p L ), vr = y/2((p R ) and vs = ($), and the matrices Z and Yp defined by (1.5) - 
(1.7) as follows: 

m{ = -^=y{v L Z , m f R = -^=y R v R Z , M F = y f s v s Y F . (3.1) 
The evolution of the Yukawa coupling constants are given by 

J t HvlZ) = ^(T/ - G{ + Hi) , (3.2) 

it HyfRZ) = vh {T «~ GfR+ hSr) ' (3 - 3) 
it HyfsYp) = itr^ ~° fs+ H * s) ' (3 - 4) 

where T* , G* and H* denote contributions from fermion-loop corrections, vertex corrections 
due to the gauge bosons and vertex corrections due to the Higgs boson, respectively. 

What is of great interest to us is to see whether the evolutions can explain the value 
R(mz) — 3 or not, i.e., our interest exists not in the hierarchy among m e , m M and m T , but in 
the hierarchy among up-quark, down-quark, charged lepton and neutrino sectors. Therefore, 
we neglect the scale-dependency of the matrix Z, because we can regard the value of z s as 
z 3 ~ 1 from Eq. (1.8). We also neglect the scale-dependency of the matrix Y F because the 
matrices Yp are expressed as 

Y F = diag(l,l,l + 36/) , (3.5) 



on the basis on which the matrix Yp is diagonal, and we find that the forms Ye =diag(l, 1, 1) 
and Yjj =diag(l,l,0) are scale-invariant and Yd =diag(l, 1, 1 — 3e tl3d ) is almost scale- 
invariant. For convenience, we approximately still use the evolution equation of y,sYp at 
/i < As and that of drZ at \i < Ar. Then, the ratio -R(/i) defined by Eq. (1.10) can be 
expressed in terms of the Yukawa coupling constants i/l, x/r and y$ as follow: 

The evolution of the ratio -R(/i) is approximately given by 



^-kiR(n) ~ 

dt vw 16tt 



\nR(tJ l )~- T7 ^(G-H) , (3.7) 



where 



H={Ht + Hl-H^)-{Hl + H%-H%) . 

(3.8) 

The G- and if -terms are given in Table 0. Since in the present model, \y L \ 2 — \yi\ 2 , |y£| 2 — 
|y^| 2 , and so on, differently from other models where \ypjyW — mf,/m t and lyl/y^l — Tn u /m T , 
we can neglect the H^- and if^-terms in Eq. (3.8). When we also neglect the i^-terms, the 
ratio i?(/i) is approximately evaluated as follows: 

/ u ,. \ 7 / 10b l 



R{m z ) 



(l + ^a 3 (m z )\n-^) 1 + °^-d \{m z ) In JL) , (3.9) 

V 2vr m z J \ 2tt m z J 



for the regions I, II and III, respectively. By using (3.9) - (3.11), we can obtain the energy 
scale \i = Ayu at which the ratio R(fi) takes R(Ayu) — 1- 
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In Fig. 2, we illustrate the behavior of Ayu for a given value of Ar. For reference, we also 
illustrate the behavior of A^°, at which a{ takes af 1 (AJ°) = 0. The value of Ayu must 
be lower than the value of A^°. Therefore, as seen in Fig. 2, if we adhere to the constraint 
&r(Ar) = cxl{Ar), we must abandon a model with a higher k value (k > 10 2 ). Only a model 
with k ~ 10 is acceptable. However, if we admit a strong coupling of the right-handed weak 
bosons at /x = Ar, for example, a^(A^) > l/An, a model with a higher k value also becomes 
acceptable. 

Of course, from a similar study, we can find that the evolution of R u /d{n) = 
{ylyjt/y^/iVLyR/ys) s ^ m keeps R u /d(mz) — 1- Therefore, the parametrization (m K,/\) u = 
(m K,/\)d in Ref. @ is justified. 



IV. EVOLUTION OF THE PATI-SALAM COLOR 

In order to avoid the burst of the U(l) gauge coupling constant, we consider that the 
U(l)y xSU(3) c symmetries are embedded into the Pati-Salam SU(4) symmetry |T3] above 



(j, = As- In other words, the SU(4)ps' gauge symmetry is broken into SU(3) c xU(l)y at 
H = As- Indeed, the structures of the heavy fermion mass matrices Mr are flavor-dependent. 
The fermions / and F belong to f L = (2, 1, 4), f R = (1, 2, 4), F L = (1, 1, 4) and F R = (1, 1, 4) 
of SU(2) x SU(2) R xSU(4) P5 at fi > A s . 

In the region III (As < \i < Ax), «l(a*) and olr{\i) are evolved with the coefficients b 1 ^ 1 
and b R n given in Table 0, but azi(fi) and a 3 (fi) are replaced with a 4 (/i) which is evolved 
with the coefficient 

b{ n = 20/3 , (4.1) 
where the boundary condition at [i = As is 

ai(A s ) = a 3 (A 5 ) = a 4 (A s ) . (4.2) 
Since oti(Ag) and a 3 (A s ) are given by 
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«r 1 (A 5 ) = \ 

a 3 \A s ) = a^(A L ) + | In + ^ , (4.4) 

Z7T Z7T A# 

respectively, the values of Ar and As are fixed at 

A^ = 5.46 x 10 12 GeV , A s = 2.73 x 10 14 GeV , (4.5) 
under the conditions (2.6) and (4.2). The unification scale A GUT is also fixed at 

A GUT = 5.84 x 10 17 GeV , (4.6) 

by the condition 

ocl(^gut) = at R (A GUT ) = a 4 (A GUT ) , (4.7) 

for example, for the embedding into SO(10) In Fig. 3, we illustrate the behaviors 

of a~ 1 (/i). Roughly speaking, the value Ar ; ~ 10 12 is favorable to the scenario B for the 
neutrino mass generation. 

On the other hand, the evolution of R(fi) defined by (3.6) is almost constant at > As, 
i.e., R(As) — R(Ayu), because there is no difference between quarks and leptons in the 
region III (As < n < Ayu)- Therefore, we obtain 

R(m z )/R(A YU ) ~ R(m z )/R(A S ) = 2.3 , (4.8) 

and we fail to obtain our desirable relation R(mz) / R(Ayu) — 3. If we adhere the unification 
of the gauge symmetries SU(3) c xSU(2) i xSU(2) R xU(l)y into a Pati-Salam type unification 
G, we must abandon the idea that the discrepancy R(mz) — 3 between quarks and leptons in 
the model given in Ref. || comes from difference of evolutions between quarks and leptons, 
or we must consider that the magnitudes of the Yukawa coupling constants are different 
between quarks and leptons from the beginning at fi = Aqut- 
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^ (Ai) + 27 ln A7 "2 



a- L \A L )+ b flu^ 
Zir A l 



2vr A R 



(4.3) 



V. CONCLUDING REMARKS 



In conclusion, we have investigated the evolution of the universal seesaw mass matrix 
model under the gauge symmetries SU(3) c xSU(2) i xSU(2) jR xU(l)y. The symmetries can 
be embedded into a unification symmetry of the Pati-Salam type at Aqut = 5.84 x 10 17 
GeV. Then, the value A R = 5.46 x 10 12 GeV is favorable to the scenario B for the neutrino 
mass generation. However, we cannot explain the discrepancy R(rriz) — 3 between quarks 
and leptons by the evolution of R(f-i) starting from R(Agut) = 1- 

On the other hand, if we abandon the grand unification scenario, the model has a possi- 
bility that the value R{mz) — 3 can be understood by the evolution of the Yukawa coupling 
constants. As seen in Fig. 2, we require cxl{A r ) = a R (Aji), the value A^ for the case which 
gives R{mz) — 3 must be A R < 10 4 GeV. If we accept a model with a strong SU(2) R force 
at /i = A R , for example, a R (A R ) = 1/47T, the region A R < 10 18 GeV also becomes allowed. 
We consider that the model with k ~ 10 is likely. Although this case rules out the scenario 
B for neutrinos, phenomenologically we can expect an abundance of new physics effects ||, 
t' production, FCNC effects, and so on, in the near future colliders. 

Our numerical results have been obtained by the constraint A R /As = 0.02 which has 
come from the observed ratio oim c /m t under the special model with (1.5) - (1.7). Since the 
ratio A R / As is fixed by the ratio m c /m t (or m^/mt) as far as a universal seesaw model with 
detMfy = are concerned, the value of the ratio A R / As is, in general, of the order of 10~ 2 . 
Therefore, our conclusions will be unchanged as far as the orders are concerned. 

In the present paper, we have not discussed a SUSY version of the present model, al- 
though the case is attractive from the point of view of the grand unification. In such a 
SUSY version, since the coefficient 803 in G-terms in (3.8) [also in Table V] is changed for 
(16/3)«3, the case push the energy scale Ayu to an unlikely ultra-high energy scale (> 10 23 
GeV). If we want to adopt a SUSY version of the present model, we must abandon the idea 
of the unification of the Yukawa coupling constants. 
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TABLES 

TABLE I. Fermion mass generation scenario 



Energy scale 


d- Sz e-sectors 


u-sector (i / 3) 


At n = As ~ XrriQ 


m(F L ,F R ) ~ A s 


m(U Li ,U Ri ) ~ A s 


At fx = A R ~ nniQ 




m(uR3, U L3 ) ~ A R 


At \x = A^ ~ rriQ 


It £ \ A L A R 
As 


m(u L3 ,U R3 ) ~ A L 

, x AlA^ 
m{u Li ,u Ri ) ~ 

A,<? 



TABLE II. Neutrino mass generation scenarios: Af± = (A^ ± N R )/\/2 



Energy scale 


Scenario A 


Scenario B 


At n = A uS 


m(N L ,N R ) ~ A uS 




At n = A s 




m(N + ,N c + )~A s 


At n = A R 




m(i/R,N-) ~ A/? 


At n = A L 


f ^ A L A R 

m{v L ,v R )~ 

Aus 


A 2 
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TABLE III. Quantum numbers of the fermions / and F and Higgs scalars <Pl, 4>R an d 3>- 
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TABLE IV. Coefficients in the evolution equations of gauge coupling constants. 





A L < n < A R 


A R <fi<A s 


A s < /i < Ax 


SU(3) C 


bi = 7 


bi 1 = 19/3 


= 3 


SU(2) L 


b{ = 19/6 


b 1 / = 19/6 


&P = 19/6 


SU(2) R 


b[ = -41/10 


b» = 19/6 


b R n = 19/6 


U(l)y 


b{ 1 = -43/6 


b[ n = -41/2 
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TABLE V. G- and if-terms in the evolution equation (3.8). i?£, Hf, and H s are given by the 
replacements u — > e and d — > f in iTjJ, and i?<j, respectively. 





I (Al < A* < a r) 


// (A fl < // < As) 


III {A s < ii < Ayu) 


G = 


47r (8a 3 - 


An (8a 3 - 2ai) 




1 (lyEl 2 - |y!l 2 ) 


H U R = 





1 (\V U R \ 


2 - K\ 2 ) 


H u s = 





m\ 2 
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FIGURES 



FIG. 1. Behaviors of a 3 (fi) (dotted chain line) , a L (fi) (dotted line), a R (/j,) (broken line) and 
a{ l (ii) [a / 1 _1 (^)] (solid line) in the case with the input values Ar = 10 5 GeV and a#(A#) = 1/Att. 
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FIG. 2. Behavior of Ayu versus Ar. The bold and thin solid lines denote the cases of the 
input values or(Ar) = ai(Aji) and au(Aji) = l/4ir, respectively. For reference, the behaviors of 
Af° for the cases of the input values or(A/j) = ai(Aji) (bold broken line) and an(Aji) = 1/Att 
(thin broken line) are illustrated. The physical value of Ayu must be Ayu < Af. 
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FIG. 3. Behaviors of qJ 1 ( / u) [aj^/i) for A5 < fi < Aqut] (dotted chain line), a^ 1 ( / u) 
[= a i? 1 ( / u) for n > Ar] (dotted line) and a^ l (n) [a' 1 _1 (/i) for A < \i < Ar] (solid line) in the 
case of the Pati-Salam type unification, where Ar = 5.46 x 10 12 GeV, A s = 2.37 x 10 14 GeV and 
Aqut = 5.84 x 10 17 GeV. 
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